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Abstract 

An SU(3) chiral Lagrangian for the lightest decuplet of baryons is constructed on 
a discrete lattice of spacetime points, and is added to an existing lattice Lagrangian 
for the lightest octets of mesons and baryons. A nonzero lattice spacing renders all 
loop integrations finite, and the continuum limit of any physical observable is iden- 
tical to the result obtained from dimensional regularization. Chiral symmetry and 
gauge invariance are preserved even at nonzero lattice spacing. Specific calculations 
discussed here include the non-renormalization of a conserved vector current, the 
magnetic moments of octet baryons, and the irN and KN sigma terms that relate 
to the nucleon's strangeness content. The quantitative difference between physics 
at a nonzero lattice spacing and physics in the continuum limit is easily computed, 
and it represents an expectation for the size of discretization errors in corresponding 
lattice QCD simulations. 
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I. INTRODUCTION AND DISCUSSION 

Chiral perturbation theory (ChPT)[l|, [J is a low momentum effective field theory 
for QCD written as an expansion in small momenta and quark masses, and it has become 
an invaluable tool for subatomic physics. With only the lightest octets of pseudoscalar 
mesons and spin-1/2 baryons, ChPT is order-by-order renormalizable and physical results 
are independent of whichever regularization prescription is chosen. 

The addition of the lightest decuplet of spin-3/2 baryons introduces a new physical scale, 
the mass difference between decuplet and octet baryons, which does not vanish in the chiral 
limit (i.e. when quark masses vanish). However, this mass difference is similar in size to 
the pseudoscalar meson masses, so the strict chiral expansion can be generalized to a "small 
scale expansion" where power counting is now meaningful even if the decuplet baryons are 
present. 0,|] 

Although physical results are independent of regularization scheme, sometimes one 
scheme is preferable over another such as when issues other than experimentally-observable 
quantities are of interest. For example, the convergence of the ChPT expression for an 
observable is often evaluated by comparing the relative sizes of the contributions that oc- 
cur at successive orders in the expansion, even though these individual contributions are 
not physically observable. The work of Refs. |5|, |6| discusses this in some detail, and goes 
one step further by associating the renormalization scale with the physical scale of baryon 
substructure. 

For another example of important non-observables, consider lattice QCD. Numerical 
simulations must be performed at nonzero lattice spacing, and computed results therefore 
differ from the desired continuum values. One might expect that the effective theory for 
a discretization of QCD is a discretization of ChPT, i.e. the most general effective theory, 
written in terms of hadronic degrees of freedom, that respects chiral symmetry and the 
other symmetries of discretized QCD and that exists in the same discretized spacetime. Our 
previous work, Ref. 0, provides one particular discretization of ChPT along with explicit 
calculations of differences between results in this theory and results in the continuum. 

The present work is mainly intended to exemplify the technique of lattice regularization 
(with explicit decuplet fields), and to note some of its features. Within a chosen lattice 
ChPT, it is a simple matter to rigorously determine extrapolations in both quark mass and 



3 



lattice spacing. Although less familiar than dimensional regularization, lattice regularization 
has the advantage that loops integrals can easily be performed numerically, since it is a 4- 
dimensional theory. The size of lattice spacing artifacts depends on the particular observable 
that is studied; at lattice spacings typical of lattice QCD simulations, we find that our lattice 
ChPT Lagrangian leads to noticeable discretization effects for the magnetic moments of octet 
baryons but negligible discretization effects for the sigma terms. 

It is important to interpret our numerical results appropriately. If one is interested in 
using a lattice chiral Lagrangian to rigorously determine the lattice spacing effects that are 
present in a specific lattice QCD simulation, then it will be important to use the chiral La- 
grangian that properly corresponds to whatever lattice QCD action was employed. Rupak 
and Shoresh|| provide a list of five 0(a) terms that can appear in the meson chiral La- 
grangian appropriate for a Wilson-type lattice QCD action, (a denotes the lattice spacing.) 
Each of these terms is multiplied by its own parameter which depends on the particular lat- 
tice QCD action of interest. For example, all five of the coefficients can be made to vanish 
by using an "improved" lattice QCD action. A similar study could be performed for the 
baryon Lagrangian. Our present work does not refer to any specific lattice QCD data, and 
this leaves the parameters at 0(a) and beyond undetermined. These parameters get speci- 
fied implicitly by our choice of a minimal Lagrangian where the only lattice spacing effects 
are contained within simple covariant derivatives and field representations. The numerical 
results for magnetic moments and sigma terms presented in this work are therefore simply 
examples of lattice spacing effects that arise from a Lagrangian which has not been specially 
improved in any way. The main point of our present work is to demonstrate the use of 
lattice regularization for chiral Lagrangian calculations, but we wish to emphasize that this 
regularization technique can certainly be applied to extensions of our minimal Lagrangian 
if one wishes to fix the lattice spacing parameters to a particular lattice QCD action. 

In Sec. II of the present work, the Lagrangian of Ref. is extended to include the 



decuplet of spin-3/2 baryons. Section [QJ contains a discussion of the electromagnetic vertex: 
the non-renormalization of this vertex at vanishing momentum transfer is shown analytically, 
and the octet baryon magnetic moments are calculated as a function of lattice spacing. In 
Sec. [IV], the KN sigma terms of the nucleon's scalar vertex are determined as a function 
of lattice spacing, with the necessary chiral counterterms determined from known baryon 
masses and the known 7rN sigma term. The running of the irN sigma term from q 2 = to 
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the Cheng-Dashen point is also determined as a function of lattice spacing. 



II. A LATTICE CHIRAL LAGRANGIAN 

The standard SU(3) chiral Lagrangian containing pseudoscalar mesons (M), spin-1/2 
baryons (B) and spin-3/2 baryons (T) is 

£ = Cm + £mb + £mt + £mbt- (-0 

To expand in powers of external momenta, of meson masses and of the T-B mass difference 
(relative to the larger scales: baryon masses and the chiral symmetry breaking scale A x ~ 
AtcF-k), it is standard practice to write the Lagrangian in terms of heavy baryon fields, B v (x) 
and T Vfl (x), instead of the relativistic fields, B(x) and T^x). The transformation for the 
spin-1/2 field is 

B v (x) = exp(im RB v ■ + f)B(x), (2) 

with mHB chosen near the average octet baryon mass. A similar transformation is used to 
define the decuplet field ||. The expanded Lagrangian becomes 

r - 4- r<® 4- 4- /" (2) 4- /" (3) 4- 4- 4- ("Vl 

where the omitted terms only contribute to octet baryon properties beyond third order. In 
Euclidean spacetime, 

(4) 



4? = ^Tr^VrrfVpU-rfU-xU^ 



£mb = ( m o -m H B)Tr [B V B V ] , (5) 

(6) 



41b = E Tr (B vVlx D, L B v ) + X?Tr (B v S„{u^, B v }) + ^Tr B v 
£mb = ^Tr (B v {v -DvD- D 2 )B^ - b v Ti (b v { X +, B v }) - b^Tr (b v [ X+ , B v 



2m 

ljJL D 



-6 Tr (B V B V ) Tr + £ Tr ^[^, S„]{£i££t + C^C, 



4mo 



+ ETr fe^, + ^if^, 5,]) + • • • , (7) 

^MB — ■ • • > (8) 



ijk 



(9) 



4& = -EEfevvr*? -wEW2S')-aee^. 

A»f ijfe V i ' V ijk 

4& T = % E E ^ (t v T u 7 b v + BX mT H lm ) > ( 10 ) 
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A* ijklm 
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where = -75 J2u a ^vV u is the Pauli-Lubanski spin vector, A is the decuplet-octet mass 
difference, and we choose = (0, 0, 0, 1) so that the covariant derivatives of Ref. |7| can be 
retained. The field strength tensors Fk,(x) and F^ v {x) correspond to external spin-1 fields. 
The superscripted indices in Eqs. (0) and (110) are flavour indices: 



(Bf 
Bf 
\Bf 



/tt 11 


TT 12 
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7T 21 


TT 22 


71 
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TT 32 
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Bf 


Bf 
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23 



33 
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32 



S 33 
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/ J_v0 1 J_a 
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111 
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Vfl 

rpll3 

Vfl 



3,3 / 

D / 
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T^123 
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V2' 



K° 



\ 



(11) 



V 



V 

rp222 
Vfl 

rp223 
Vfl 

m233 

Vfl 

T^333 
vfj, I 



Vfl 



V2 « + V6 v 

wo 

' — 'v 

-^A+ -^A° 

V)j. ^/3 vfi 



-fa ) 

Pv 

A 



(12) 



'Vfl 



1 T*+ 



J_y*0 

S^ufi, 



\ 



1 v*— 
73 V M 



V3^M 



(13) 



where Tf^ is understood to be completely symmetric in i,j,k. In the Lagrangian, the 
pseudoscalar mesons are represented nonlinearly, 

U(x) = £ 2 {x) = exp{-i\ a 7c a {x)/F), (14) 

where A a is a Gell-Mann matrix, and the current quark mass matrix, A4, enters via 

X = 2BM, (15) 



X+ 



(16) 



So far the definitions of this section have been general, and have not assumed a spacetime 
lattice at all. We now consider the derivative structures in the Lagrangian, and it is here that 
the expressions become lattice dependent. As in Ref. @, we choose an isotropic hypercubic 
lattice. Nearest neighbour sites are separated by a distance a, and a M will be used to denote 
a vector of length a in the positive /1 direction. 

First, we recall the non-decuplet definitions from Ref. 0. Denoting the external spin-1 
fields corresponding to Fz,(x) and F^ v (x) by L^{x) and R^x) respectively, the derivatives 



that appear in the lowest-order meson Lagrangian, Eq. 

1 

a 



V<%(x) 



R IJi {x)U(x + afjL^x) 



are taken to be 



U(x) 



(17) 
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and the meson derivative needed for the baryon Lagrangian is 



where 



V&U(x) = — [R,(x)U(x + a M )Lt (x) - Rl(x)U(x - a^L^x - a 
The covariant derivative for the octet baryon field is 



aD i B v (x) = B v {x) — r 4 ' (x — a^B v {x — a^)Y\{x — a 4 ), 

x,y=l.r 



aDjB v (x) — — ^2 [rf(x)B v (x + aj)rJ\x) — Tf\x — aj)B v (x — aj)r^(x — aj) 

^ X,Y=L,R 



(18) 



(19) 



(20) 



(21) 



in the temporal and spatial directions respectively, where 



C,(x)L ll (x)^(x + dp), for X 



(22) 



^(x)R ll (x)^(x + a lt ), ioiX = R. 

Not defined explicitly in Ref. JF| but required for the present work are the field strength 
tensors, which we discretize as follows so that the chiral transformation properties and the 
lattice symmetries are respected: 

Aia 2 F*(x) = A-X^x)X v (x + a^)Xl(x + a u )Xi(x) 

-X v (x)Xj t (x - a fl + a v )Xl(x - a v )X^(x - a M ) 

-Xp(x - a^)Xl(x - dp - d v )X^(x - a M - a v )X v {x - a v ) 

-Xl(x - d v )X^(x - a v )X v [x + a M - a^Xfa), (23) 

for X = L and R. Finally, consider the covariant derivative of the decuplet field. Since we 
have chosen = (0, 0, 0, 1), only the temporal derivative is needed and we define it to be 

dD 4 Tti k (x) = THHx) - [Fii(x - a 4 ) + r* (x - a 4 )] Tg k (x - a 4 ) 



■ Vfl 



1 

J2 - 04) + r 4w(* - «4)] - 04). 



(24) 



The chiral invariance of this Lagrangian can readily be verified by adding the chiral trans 
formation of 
given in Ref. 



formation of the decuplet field, T l J k — > ouOj m Ok n T^ n , to the set of chiral transformations 
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When performing calculations, it is important to remember that the decuplet propagator 
must always be accompanied by the projector which eliminates spurious spin- 1/2 compo- 
nents from the propagating field. For the present Lagangian, the product of propagator 
and projector is 

Ttt sin(ap 4 ) — i[aA + 2 sin 2 (ap 4 /2)] 

III. GAUGE INVARIANCE AND OCTET BARYON MAGNETIC MOMENTS 

To calculate the electromagnetic form factors of an octet baryon, one simply identifies 
the photon field, A^x), with the external spin-1 fields of Sec. II, 

L^(x) = R^{x) = exp [-iaeQA^x)] , (26) 

where Q = diag(2/3, —1/3, —1/3). 

From Eqs. (j3UT0l), it is easy to see that the leading order contribution to an octet baryon's 
electromagnetic vertex is contained within C^ B and the next-to-leading order contribution 
within £^b- The resulting matrix element is 

B(p')\j™\B{p)) = zeS 1 ,(A; / )[^gi ? ^l-2sin(aA; 4 )-2sin 2 ^^ 

LO 



+ Y1 S » S »1 sin ( a ^) f 1 - \ s in ( a ^) - sin 2 i^f] 

(XfTlQ u \ Z \ Z / 



where q = p' — p is the momentum transfer, Qb denotes the electric charge of the baryon, 
and fi B ° are the leading order expressions for the magnetic moments: 

(28) 
(29) 
(30) 
(31) 
(32) 

Notice from Eq. (27) that the term containing Qb is purely temporal and the term is 
purely spatial with our chosen frame, since — (0, 0, 0, 1) =>■ 5" • v — 0. 

Returning to the full Lagrangian of Eq. (|3|), one finds the corrections to the matrix element 
of Eq. (p7|) that are shown diagrammatically in Fig. |l| As will now be shown, evaluation of 
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FIG. 1: One-loop contributions to an octet baryon's electromagnetic vertex. Dashed, solid, double 
and wavy lines represent mesons, octet baryons, decuplet baryons and photons respectively. 5Z 
denotes the contribution to wave function renormalization that arises from the diagrams in Fig. ^. 

these diagrams leads to the renormalization of /i# , an d to the non- renormalization of Qb 
as required by gauge invariance. 

A. The non-renormalization of electric charge 

To verify the non-renormalization of Qb-, it is sufficient to work at vanishing momentum 
transfer. Also, the momentum of each external baryon is simply (0, 0, 0, itlb) plus corrections 
which are of negligibly high order in the chiral expansion. For definiteness we will discuss the 
proton's charge in this subsection; the extension to other octet baryons is straightforward. 

In this limit, the contribution of the diagram in Fig. |I|(a) to the matrix element in Eq. fl27|) 
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(a) 



(b) 



(c) 



i 

\ 



_fc 



FIG. 2: One-loop contributions to the wave function renormalization of an octet baryon. Dashed, 
solid and double lines represent mesons, octet baryons and decuplet baryons respectively. 



with a 7r + in the loop is 



M 



(a) 



t/o d 4 g f —iev 



a 



-./a (2tt) 4 V 2F 2 y V 2 Ea[1 -cos(ag A )] + a%2 
-\ ,'„/■-, /••■••• ' i </ 

exp 



-ia^ev^ f°° f K / a d 4 g 



4F 2 



o 



yr/a (27r) 



-2 (4 + ^ - ^ cos(ag A ) 



-iev 



»WAa 2 x 2 ) 



Aa 2 F 2 

where xm is related to the meson mass[|7| 

2 . /ax M \ 

m M = -arcsmh , 

a \ 2 / 

and Wi{e 2 ) is the integral transform of the fourth power of a Bessel function 0, 



W n (e 2 ) = / dx/ n (x)exp 
Jo 



-x I n + — 



(33) 



(34) 



(35) 



Summing over all internal mesons gives 



(a) 
Af 



—lev^ 
Aa?F 2 
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^) W^xl) + 2 [l - ^) iy 4 (a 2 * 2 



(36) 



Notice that this contribution is quadratically divergent as a — > 0. Dimensional regularization 
does not show power divergences, and Fig. |I|(a) vanishes exactly in that scheme. 

Figure |l](b) is nonzero only if the loop meson is 7r =1= or K . For the pion loop, the lattice 
regularized expression is 



M 



ip) 



iea 2 v^ f n / a d A q sin 2 (ag 4 



4F 2 J-v/a (2vr) 4 [4 + a 2 x\ /2 - £ A cos(aq x )Y 



4a 2 F 2 
-ieVfj, 
4a 2 F 2 



dzl*{z)l' {z)&cp 



-z 4 + 



a 2 x 2 ' 



a 2 x 2 



22- 



1 + — ^ W 4 (a 2 x 



(37) 
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Repeating this calculation for the kaon loop leads to 



M {b) =Y, M 



M 



M 



—lev^ 
4a 2 F 2 



3 + 

4 I 8 



W 4 (a 2 xi) -21 + 



WJax 



2^,2 \ 
Kj 



• (38) 



The only other contribution to the renormalization of Qb that is independent of the axial 
couplings {V and T) and the decuplet (with coupling C) comes from the tadpole contribution 
to wavefunction renormalization. As discussed in Ref . , the lowest-order octet baryon two- 
point function is 



bb = m KB - m 



sin^aq^) — 2% sin' 



\ 2 



The one-loop correction shown in Fig. |2](a) evaluates as follows: 



pp 



—i 
a 



sin(afc 4 ) — 2i sin 2 



afc 4 



6Z, 



(a)' 



where the external proton momentum is m$v + k and 



p 3a 2 F 2 



1 - ^a 2 xlW A {a 2 xl) - ^a 2 x 2 K W A {a 2 x 2 K ) - ^a 2 x\W A {a 2 x%) 



(39) 



(40) 



(41) 



This is precisely what was required to facilitate the expected non-renormalization, since 
Eqs. (H), (H) and (P) give 



M (a) +M (b) +iev fl 5Z^ = 0. 
Turning to the contributions that depend on T> and J 7 , one finds 



M 



(c) 



M { v c) 



where 



G M (e 2 ) 



iea 2 v^ f 77 / 11 d 4 g 



(V-F) 2 G (a 2 x 2 K ), 
i(D-3^) 2 G (a 2 4) 



[cos(ag 4 ) — i sin(ag 4 )] J2l=i sin 2 (agfc) 



(42) 

(43) 
(44) 
(45) 



AF 2 J-n/a (2vr) 4 [4 + e 2 /2 - Ea cos(ag A )] [sin(ag 4 ) - iaM - 2z sin 2 (ag 4 /2)]< 



(46) 



The contribution of Fig. |T](d) with a charged pion in the loop is 



m«! = — ( v + rfj 



n/a d 4 q 



Vn sin(ag 4 ) J2 P S p sm(aq p ) £ ff S a sin(ag CT ) 



iea 2 v 
2F 2 
d 



^(V + F) 2 x 

Tr/a d 4 g 



/a (2tt) 4 [sin(ag 4 ) - 2i sin 2 (ag 4 /2)] [4 + a 2 x\ /2 - £ A cos(ag A )] 
sin(ag 4 ) ELi sin2 ( a <?fc) 



d(a 2 x 2 ) J--K/a (2tt) 4 [sin(ag 4 ) - 2% sin 2 (ag 4 /2)] [4 + a 2 x\ /2 - £a cos(ag A )] 
(P + ^) 2 Go(a 2 x 2 ), 



(47) 
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where we have used S4 = and S 2 = S* 2 = Sf = —1/4 as well as integration by parts. 
Similarly, 

Mf + = \iV 2 + ^ 2 )G,(a 2 x 2 K ). (48) 

Next, consider the wave function renormalization of Fig. 0(b). For a tc° loop, the contri- 
bution to the two-point function is 

rmm(-!T /a,Tr /a— k) (\^q 



X 



— min(7r/a,7r/a+fc) (271")^ 

E p S p sin(ag p ) J2a S a sm(aq a ) 



[sin(ag 4 + a/c 4 ) — 2i sin ((ag 4 + a/c 4 )/2)][4 + a 2 xl/2 — J2\ cos(ag A )] 
sinfa/ci) — 2i sin 2 (a&; 4 ) — a5X 



SZ® L 



/a 



(49) 



where the external baryon momentum is m^v + A; and <5X produces a mass renormaliza- 
tion. The limits of integration have been chosen to ensure that the momentum of each 
internal propagator remains within the lattice's Brillouin zone over the entire domain of q 
integration. As it happens, the ^-dependences in the limits of integration do not affect 
wave function renormalization, and the result is 

iev.SZ® =~{V + F?G Q {a 2 xl), (50) 

which exactly cancels the 7r° loop from Eq. ([43]) . In the same way, replacing the 7r° in SZ 
by each of the other mesons serves to exactly cancel Eqs. (|44]) , (|45|) , (f47|) and ([48l). 

Finally, the decuplet contributions of Figs. |l](e), |l](f) and 0(c) are found to cancel in 
essentially that same manner as did the octet contributions: 

Mil = ^ 2 G A (a 2 * 2 ), (51) 

*e5Zl% = -jjc 2 G A (a 2 * 2 ), (52) 

M$ = -ie5z£ = ^C 2 G A (a 2 * 2 ), (53) 

M%\ = -ie6Z$ = lc 2 G A (a 2 x 2 K ), (54) 

9 

M { H = -^C 2 G A (a 2 * 2 ), (55) 

M ( il = -iebZ% = l -C 2 G A {a 2 x 2 K ). (56) 

Thus, all corrections to the proton's electric charge at first loop order in the chiral expansion 
sum to zero, and the non-renormalization of electric charge is confirmed. 
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B. The octet baryon magnetic moments 

The one-loop corrections to octet baryon magnetic moments come from Figs. [l](d) and 
|l](f). These diagrams were discussed in the previous subsection, but only for vanishing 
momentum transfer which is not sufficient to obtain magnetic moments. 

It is convenient to choose the Breit frame where the incoming baryon has momentum 
iTL v — q/2 and the outgoing baryon has momentum m v + q/2, and to choose the integration 
momentum to be the internal baryon's momentum. This choice displays symmetries in the 
integrand that help to simplify the calculations. Notice that higher order corrections to the 
external baryon momenta have already been omitted, and that mo is equal to the physical 
mass in these loop diagrams since the difference is of higher chiral order. 

The evaluation of Figs. |l](d) and |l|(f ) for arbitrary momentum transfer q gives 

-2e 

am^v 

where 8F2 is the loop correction to the Pauli form factor. The magnetic moments are 



+ = — £0,., S v \ sin 5F 2 (q), (57) 



IL B = v h B° + 6F?*{0) B + SFt c (0) _ , (58) 



B 



?dcc 



B 



where the last two terms are from Figs. |l](d) and 0(f) respectively. Both of these terms rely 
on a single integral with u 7^ 4, v 7^ 4 and fi 7^ u, namely 

-ia 3 m N r w / a d 4 g sin 2 (ag At ) cos(aq u ) 



2 -ia"m N r 



-Tr/a (2vr) 4 [4 + e 2/2-£ A cos(ag A )] 2 sin(ag 4 ) - iaM - 2i sin 2 (aq A /2) 

(59) 

but with one key difference. For the decuplet diagram it is that enters, and the mass 

splitting A ensures that the integrand has no singularities on the domain of integration. In 

that convenient form for numerical evaluation is 

2 a 3 m N f 7T / a d 4 q cos(ag Ai ) cos(ag^)[l + aA — cos(ag4)] 

M e ) = 4F 2 j_ n/a (2vr) 4 [4 + e 2 /2 - Ex cos(ag A )] [(1 + aA)(l - cos(ag 4 )) + a 2 A 2 /2] ' 

(60) 

For the octet diagram there is a singularity, and we must integrate around it according to 
the usual "+ie" prescription for field theory. Some details of this procedure in the context 
of lattice regularization can be found in the appendices of Ref. M, so here we simply provide 
our final expression: 
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TABLE I: Coefficients that appear in the octet baryon magnetic moments. 
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^A 1 d 3 g cos^g^) cos(aq v 



-n/ a (2tt) 3 [3 + e 2/ 2 - £j* =1 cos(ag* 

This too is easily evaluated numerically. 

Our final expression for the octet baryon magnetic moments is 



(61) 



Vb = cffir, + 4> F + c%H (a 2 x;) + c%H (a 2 x 2 K ) + ct*H A (a 2 x 2 ) + c%H A (a 2 x 2 K ), (62) 

where the coefficients are listed in Table |. Hm diverges as a — > 0, but the offending terms 
can be absorbed into renormalized definitions of fin and /j,p. The result is identical to that 
obtained from dimensional regularization. However lattice regularization also allows us to 
compute at nonzero a, and it is interesting to consider lattice spacings that are typical of 
lattice QCD simulations. 

Using the experimental values of \x v and fi n to fix the parameters \in and fip, the other 
seven magnetic moments become predictions of the theory, and are plotted as a function 
of the lattice cutoff n/a in Fig. |3]. The plot assumes standard experimental values for the 
coefficients appearing within loops: V = 0.75, T = 0.50, and C = 1.5. These values could 
be varied within experimental uncertainties, but such details do not significantly affect our 
present interest: the size of discretization effects. 

As seen in Fig. |^, each magnetic moment smoothly approaches the corresponding dimen- 
sional regularized result at n/a — ► 00. These limiting values are given in Table [TT]. The 
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FIG. 3: Octet baryon magnetic moments (in units of //jy) as functions of lattice spacing. La- 
grangian parameters are fixed by requiring the proton and neutron magnetic moments to equal 
their experimental values at all lattice spacings. 



TABLE II: A comparison of the magnetic moments at ir/a = 6.0 GeV, corresponding to a = 0.10 
fm, and their values in the continuum limit. 



B 


£+ 


S° 


S - 


CO 




A 


AS 


H B (a = 0) 


1.64 


0.12 


-1.40 


-0.14 


-0.98 


-0.12 


1.11 


= 0.1fm)///B(a = 0) 


1.11 


1.94 


0.97 


2.58 


1.01 


1.95 


1.05 



agreement with experiment is not particularly impressive at this chiral order, as has been 



known for some time||. The situation is dramatically improved at next chiral order [M], [I I 
or by the methods of Refs. 0, |6|. At present, we focus on the discretization effects rather 
than a precise comparison to experiment. 

For ir/a = 6 GeV, Table [Tl] shows that the relative sizes of discretization effects vary 
from a few percent to a factor of 2 or more, depending on which magnetic moment is 
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chosen. The large variation is somewhat misleading: the absolute discretizations are quite 
comparable for all magnetic moments, as is evident from Fig. |3|. However, even 0(10% — > 
30%) discretization uncertainties are significant in this context, since chiral corrections are 
typically of this order too. If ChPT is being employed as a way to determine chiral effects 
in lattice QCD, then these discretization effects must be considered. Notice that it /a = 6 
GeV corresponds to a lattice spacing of 0.1 fm, which is typical of modern lattice QCD 
simulations. 



Sigma terms are the scalar form factors of a baryon multiplied by the quark mass, 



where rh = (m u + md) /2. The sigma terms vanish in the chiral limit and are therefore useful 
in discussions of chiral symmetry and its breaking. Furthermore, they offer a probe of the 
nucleon's strangeness content, 



IV. THE ttN AND KN SIGMA TERMS 




(63) 



(64) 



(65) 



y = 




(66) 



and this quantity continues to be of great interest to many researchers. (See Ref. || and 
references therein.) 

The lattice-regularized Lagrangian leads to 

3 




3=1 4> 



4Jv = -2xUb + b D )+xlY,Y, b t K A<f>A)-^(V 2 -2VF- 3^)^,7)) 



3=1 <t> 



+^K 2 (7r, V ), (68) 

3 2 

a p N = -24(i -M + 4EEft(W) + Tp 2 - 2 ^-^i(M) 



3=1 <t> 
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TABLE III: Coefficients that appear in the nucleon sigma terms. 



■ A <t> 
3 9 a i 


j 




1 vr (3/2)(V + F) 2 


(3/4)(V + F) 2 


(3/4)(V + F) 2 


1 K (1/6) (5£> 2 -6VT + 9F 2 ) 


(7/6)P 2 (5/2)J^ 2 


{3/2){V-T) 2 


1 77 (1/18)(V - 3T) 2 


(5/36) (P - 3.F) 2 


(5/36) (V - 3J 7 ) 2 


2 vr 3/2 


3/4 


3/4 


Z A o/Z 


o/z 


o/z 


2 7/ 3/4 


3/2 


25/36 


3 vr 2 


1 


1 


3 if 1/4 


1/3 


1/2 


3t/ 









where the integrals are 

, 1 /- d 4 # ELi sin 2 9 k 



16aF 2 J-* (27r) 4 [4 + a 2 x\ /2 - Ea cos a ] [4 + a 2 x 2 j2 - Ea cos 9 



A 



, _JL_ r &e ELi sin 2 e k 

16aF 2 (2tt)3 [3 + a 2 x 2 j2 - ELi cosfl fc ][3 + a 2 x 2 2 /2 - Eti cos 9 k ] ' 1 1 

K U .s 1 r d 4 fl (1- cos 6) 4 ) ( ] 

2{<Pu<P2) 8aF 2 i-.(27r) 4 [4 + a 2 x 2 i /2-EACOs^][4 + a 2 x 2 2 /2-EACOs^]' 1 ' 



-C 2 dW (l + aA-cos^ 4 )Efc=iSin 2 ^ 



A 



12aF 2 (2tt) 4 [4 + a 2 x 2 /2 - Ea cosfl A ] 2 [(l + aA)(l - cos 4 ) + a 2 A 2 /2] ' 

(72) 



and their coefficients are defined in Table [TLT| . 

Following Ref. 0, the octet baryon masses and ct w n(Q) are constrained to their physical 
values at all lattice spacings. These five observables are functions of five parameters: mo, 
bo, bx>, bp and T>. The other axial coupling is T = 1.267— T>. With these parameters fixed, 
predictions are obtained for the KN sigma terms as shown in Fig. f|. Notice that the KN 
sigma terms approach their continuum values very quickly: at nja = 1 GeV the difference 
is 0(20%) and at vr/a = 6 GeV the difference is 0(1%). 

The momentum dependences of the scalar form factors are parameter-free, so all lattice 
spacing effects are exclusively from loop diagrams. For definiteness, consider the running 
of the ttN sigma term to the Cheng-Dashen point, q 2 = —2m 2 in Euclidean notation. 
This is obtained from the scalar vertex with incoming momentum q = iQ/a where Q = 
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FIG. 4: Kaon-nucleon sigma terms (in units of GeV) as functions of lattice spacing. Lagrangian 
parameters are fixed by requiring the octet baryon masses and <7„-jv(0) to equal their experimental 
values at all lattice spacings. 



(0, 0, y/2ax n , 0). Numerical results can be obtained directly from the integral expression, 

a nN (-2xl) - <7 WJV (0) =4EE 4 [ k M) ' K M <!>)] > (73) 

j'=i 

with 

d 4 6> ELi [ sin2 °k cosh 2 Q fc + cos 2 6 k sinh 2 Q fc ] 



1 d 4 ( 



+ 



(2tt) 4 [(4 + a 2 a; 2 /2 - Ea cos 6 x cosh Q fc ) 2 + (ELi sin 6 k sinh Q fc ) 2 ] 
71 d 3 6> ELi [ sin2 ^ cosh 2 Q k + cos 2 6> fc sinh 2 Q k ] 



M4>) 



KM 



16aF 2 J-n (2tt) 3 [(3 + a 2 x 2 /2 - ELi cos 6 k cosh Q k f + (ELi sin ^ sinh Q k f] ' 

(74) 

i r d 4 e a-cos0 4 ) 



8aF 2 (2tt) 4 [(4 + a 2 x 2 /2 - Ea cos 6 x cosh Q k f + (ELi sin 9 k sinh Q fc ) 



C 2 r d 4 6 



ELJsin 2 9k cosh 2 Q k + cos 2 6» fc sinh 2 Q k ] 



12aF 2 7-. (2tt) 4 [(4 + a 2 x 2 /2 - Ea cos 6 x cosh Q A ) 2 + (ELi sin 6 k sinh Q fc ) 
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FIG. 5: The difference (in units of MeV) between the pion-nucleon sigma term at the Cheng-Dashen 
point and at q 2 = 0, as a function of lattice spacing. 



'1 + aA - cos6» 4 ) 



(76) 



(1 + aA)(l -cosfl 4 ) + a 2 A 2 /2' 

and the results are plotted in Fig. In this case, lattice spacing effects are at the few 
percent level for n/a = 6 GeV. 
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